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Abstract. The general features and characteristics of Kapteyn series, which 
are a special type of series involving Bessel function, are investigated. For many 
applications to physics, astrophysics, and mathematics, it is crucial to have closed- 
form expressions in order to determine their functional structure and parametric 
behavior. Closed-form expressions of Kapteyn series have mostly been limited 
to special cases, even though there are often similarities in the approaches used 
to reduce the series to analytically tractable forms. The goal of this paper is to 
review the previous work in the area and to show that Kapteyn series can be 
expressed as trigonometric or gamma function series, which can be evaluated in 
closed form for specific parameters. Two examples with a similar structure are 
given, showing the complexity of Kapteyn series. 



1. Introduction 

Kapteyn series of the first and se cond kind have arisen in a large variety of physics 
problems since their discovery by iKapteyr] ( 18931 ) . Series of the first kind are of the 
form 

oo 

Kx = Ki({a n },a,/3,c,b) = ^ a-nJan+picn + b) (1) 

n— — oo 

with a,/3,c, and b fixed but possibly complex, and {a n } is a sequence of complex 
coefficients. Series of the second kind are of the form 

K 2 = K 2 ({a n }, a, (3, 7, e, c, 6, /, g) 

OO 

= ^2 a n J an+0 (cn + b)J^ n+e (fn + g) (2) 



n— — 00 



with a, /3,j,e,c,b, f, and g fixed but possibly complex and {a n } is a sequence of 
complex coefficients. Over the years, the subject of summing Kapt eyn series has 



attracted both physics (see below) and mathematics researchers fe. g.. lErdelvil 11981 
WatsorJ . [l96(l . 



One of the more important aspects of such Kapteyn series is to provide closed 
form expressions for the series in particular problems of physical interest. This 
aspect is often crucial, as one attempts to determine the functional structure and 
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parametric behavior for the problem at hand. For exa mple, for a syste m of N equally 
charged particles uniformly spaced on a rotating ring, iBuddenl ([19241 ) used Kapteyn 
series to show how the far-field radiatio n distribution varies as N — > oo (see also 
ILerche and Tautzl . 120081 ). Earlier, Schott ( 1912 ) had discussed radiation from a single 
relativistic particle moving in a circle. Both radiation problems involved Kapteyn 
series of the second kind with a = 7 = 1; ft = e = 0; c = / and b = g = but with 
different coefficients a n for the two problems. 

To name but a few illustrations from a broad range of physical applications (see 
Tautz and Lerche . 20091 . for a review), the large varie ty of physics problem s involving 
various Kapteyn series includes: (i) Kepler's problem dDattoli et al.l . ll998l ) (ii) pulsars 



( Harrison and Tademaru , 19751 : Lerche and Tautzl. 12007); (iii) s ide-band spectra of 

" 20091); Ci v) high-intensity 
v) queueing 



tera-hertz electromagn etic waves (|Citrinl.ll999t ILerche et al . . 
Compt on scattering (lHarvev et al.T " 20091 Lerche and Tautzf " 2010l); 
theory (|Dominicil. 2007? ) ; and (vi) cosmic ray transport theories ( Tautz and Dominici . 
2Q10 HTautz and Lerche l2010l). 



In iDominicil (|2007a ) , Kapteyn series of the form 

00 

K 1 (z,t) = jy i J n (nz) 



n=l 



(3) 



were studied and a series representation was derived in powers of z. Furthermore, the 
radius of conve rgence wa s analy zed. Also, general Kapteyn series of the first kind were 
considered. In IDominicil (|2010j ). an asymptotic approximation in terms of a Kapteyn 
series was obtained for the zeros of the Hermite polynomials. 

As far as can be determined, physical applications of Kapteyn series to date seem 
to involve only structural behaviors of the form 



' a n J n+v (cn) 



(4) 



for Kapteyn series of the first kind and 

a n J n +p{cn)J-{n+e){cn) 

for Kapteyn series of the second kind. What 
broad ranges of the parameters and coefficients a 
representations of the Kapteyn series K\ and K2 are available. Such knowledge would 

then obviate having to evaluate each application of Kaptey n series de nov o. 

Efforts in this general direction have been provided by iNielsenl (|l90lL Il904 ) who 
sum med particular K apteyn series of the second type. Curiously, in respect of Nielsen's 
work, IWatsonl (|l966l) remarks "series of the type 



(5) 

would be of value is to determine 
so that rather general closed-form 



PnJi/-\-n 



+ n 



1 H+n 



+ n \ z 



have been studied in some detail by iNielsen (1901). But the only series of this type 
which have, as yet, proved to be of practical importance, are some special series with 
U = v, and w ith simple coefficients." However, Watson also goes on to say that 
Schottl (|l912h has shown that" 



71=1 



J n {nzf 



(1 - z 2 ) 



-1/2 



- 1 



but direct inspection of equation (31) from INielsenl (|190ll) shows that the formula 
ascribed to Schott was already available. Indeed many further direct summations of 
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Kapteyn series of the second kind are to be found in Nielsen's (|190ll ) work such as 
[his equation (31a)] 



2 E J ' 



{{2n + l)x]J' n [{2n + l)x] 



1 r 
2~r 



-1/2 



- 1 



(6) 



The difference in philosophy between Nielsen and Schott is that Nielsen treated 
the summation as a pure mathematics' problem requiring summation, while Schott 
worked out the physics problem of synchrotron radiation involving the series. Thus, 
physics applications of the Kapteyn series arose a decade (or so) after the series was 
originally summed in closed form. 

It would seem that to prejudge the ability to provide summations of as many 
as possible Kapteyn series of the second kind as not of practical use is basically not 
appropriate, for applications often follow much later than the basic mathematical 
results. 

For these reasons, it seems relevant to consider de novo the series K% and to 
attempt to determine procedures for evaluating such series for as broad a range of 
parameters and coefficients as possible. In Sees. [5] and G2 general methods as well as 
two specific examples for the summation of Kapteyn series of the second kind will 
be presented, respectively. Sec. 0] provides a short summary and a discussion of the 
results. 



2. Methods for Summing K 2 Series 

In this section, different approaches are investigated that have proven useful for 
summing various Kapteyn series of the second kind. 

2.1. An Integral Representation Procedure 

One of the main techniques for summ ing Kapteyn series of the second kind is (|Erdelvil . 
19811: iGradshtevn and Rvzhikl . l2000h 



2 



tt/2 



Jn(z)Jv(z) = - / d0 J M+w (22cos^)cos(/z- v)6, (7) 



which is valid for fi, v any integer values and is otherwise valid for 5R(/i + v) > — 1. 
Note here the requirement that J M and J„ have the same argument. 
Thus, when considering 

00 

K 2 = ^2 a n J an+ {3(cn + b)J 7n+e (fn + g) (8) 



n— — oo 



one restricts the evaluation using equation ([7|1 to c = f and b = g. Then, for 
3?(/i + v) > —1 one requires 5ft[(a + + f3 + e] > —1 for all integers n, with a n 7^ 0. 
If a n ^ for some integer n, then one requires 3?(a + 7) = and %t(f3 + e) > —1. 
Hence, one has to evaluate 

00 

K 2 = E a n J an+ p(cn + 6)J_ QR „ +i7l „ +e (cn + 6), (9) 

n— — oo 

where 7 = jr + iji with 7r = -«r and a = q:r + iai, and 9?(/3 + e) > —1. 
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Starting with the case /j, = — v in equation ([7]), one has for the Kapteyn series of 
the second kind, K 2 ({a n }, a, /3, —a, — j3, c, 6, c, 6), i. e., 

oo 

K 2 = ^2 a n J an+l3 (cn + b)J_( an+l3 )(cn + b) (10a) 

71 — — OO 

2 r /2 



2 pit / ^ 

= -/ a„J [2 cos 0(cn + b)] cos[26>(cm + /?)] 

Using the fact that 



(106) 



71 — — OO 



Jo(z) = — I dt/j cos(zsinV') (11) 



one gets 

r n/2 



K 2 = ^ I dO dip a„ cos[26>(cm + /3)1 cos[2 cos 6> sin ^>(cn + 6)1 

= 3 / <W ^ E a « 



(12a) 



n— — oo 

X 



^2 j {cos{26 (an + f3) - 2j cos8 sin ip(cn + b)}} . (126) 
To the extent that one can sum expressions such as 

OO 

s = E a n cos(nA + B) (13) 

n— — oo 

in closed form, if 2 can be, at worst, reduced to a double integral and, at best, can 
be evaluated in closed form. The reduction depends precisely on the functional forms 
chosen for a n and on the convergence of the terms in the series S (normally, but not 
necessarily, by taking A to be real). 

For example, if a n = (n + p)~ 1 where p is not an integer, then one can write 

<9S" 

S = -pS'cosB + — sin.B, (14) 
a A 



where 



OO , 

S'= V (15) 
^ n 2 - p 2 v ; 



n— — oo 



Because S' can be given in closed form, it is possible to sum a variety of Kapteyn 
series of the second kind with this procedure. 



Thus, for sideband spectra in the tera-hertz regime one can show (jLerche et al 

2009f) that 

-j— p = (-1) ^ csciKbjJn+biabjJn^ab) (16) 

V— 1 

for n integer and n ^ 1, with < a < 1 and < 6 < 1 together with 

E / I 2 u9 J n+,+i (a(v + §)) J n -u {a(v + i)) 
= (-1)" J sec(7r6)J n+ i +6 (a6)J„ + i_ 6 (a6) (17) 



Summing General Kapteyn Series 



5 



for n integer and n 1, with < a < 1 and < b < i. 

Similarly, for high-intensity Compton scattering, series arise of the form 

Y T-^r J n(na) (18) 
^ (a + n)i ™ v ; v ' 

with p and q positive integers, which can be reduced to analytic closed form apart from 

a single elliptic integral t hat has to be added to the rest of the closed form expressions 

( Lerche and Taut j . Boioh . 

Besides such evaluations where // = —v. the more general case with 5i(/i+^) > — 1 

needs to be considered. Writing equation (JTJ) with 

/i = an + ft (19a) 
v = — QfR« + vyin + e, (196) 

one obtains 

/! + v = i(7i + ai)n + P + e (20a) 
fj, — v = 2aR,n + i(ai — 71)71 + /3 — e. (206) 
Thus, one can write 

2 /•7T/2 OO 

K 2 = - d6 V a„J in(7l+ctl)+ ^ +e (2cos6'(cn + 6)) 

T JO 

w n=— 00 

x cos[6>(2a R n + i(a:i-7i)n + /3-6)] (21) 

with 3?(/3+e) > — 1. The cosine factors in equation (l2~Tj) converge if and only if a\ =71, 
so that 

2 /-ir/2 00 

K 2 = ~ I d9 a n hinai+p+e (2 cos 9{cn + b)) cos [9{2a R n + /3 - e] . (22) 

•'0 n=-oo 

Therefore, 



2 /*^c . 7fi/ \ 

Jj/(x) = — / df sin ( a; cosh t jcoshi/i, (23) 

for a; € R, so that c and 6 are real. But if v contains an imaginary part proportional to 
n, as in equation (|22[) . then the series in equation (|22p diverges exponentially unless 
a„ converges fast enough (e.g., a n proportional to exp[— bn 2 }). Thus, one requires 
on = 0. 

Under such conditions one can write 

'2\2 />oc /' 7r /2 



Ko = 



dt dO a„ cos[i^(2aftn + /3 — e)] 

Jo Jo — 



n— — 00 

7T 



x cosh [(/J + e)i] sin 2 cos 6»(cn + 6) coshi -- ((3 + e) . (24) 

Again, one sees that for choices of a n such that the series in equation (|24[) can be 
summed in closed form then K 2 is reduced, at worst, to a double integral and, at best, 
can be evaluated explicitly. 

All of these procedures for summing the general second-order Kapteyn series 
represented by K 2 are dependent on the integral representation from equation (JTJ) for 
J l/ (z)J fl (z) 1 valid for 5i(/i + ;/) > —1 when /i, v are not integer and otherwise generally 
valid. 

But just because there are values of jJi(/i + v) ^ — 1, for which equation ([7]) cannot 
be used, does not mean that other Kapteyn series of the K 2 form cannot be summed. 
One needs other procedures to effect the summations when 5R(/i + v) ^ —1. 
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2.2. Series Representation Procedures 



T he use of series representations to sum if 2 types of Kaptey n series was alr eady known 
tolNielsenl (Il90lh and later the same procedure was given bv lWatsonl (|l966l ). Following 



Nielsen! (|1901| ). the sense of the argument is as follows: One considers first integrals 
of the form 

tt/2 

(25) 



l v {a) — dx 
Jo 



cos 1 ' 1 x cos ax 



and, for 



{v) > 0, expresses the result as 
7r cos an/2 



T v (a) 



(26) 



2 u v T((v+l + a)/2))T((u + l-a)/2)' 

Then one considers equation ([7]) with fi = n — a and v = n + a. According to Nielsen, 
one then uses the series expansion definition of the Bessel function under the integral 
sign as 



J a (z) 



{-If 



.,2 k 



2 J ^ 2 2k kl T(a + k + l) 



(27) 



and so one integrates equation ([7]) term by term with [i and v as defined above by 
using equation (|27|) . Effectively, one trades a sum over Bessel functions of the Kapteyn 
kind for a power series. The resulting power series can often, but not universally, be 
either summed in closed form or can be evaluated for specific parameter values. In 
this way, Nielsen argued that 

n!T(n+ \){2x) 2n 



sini^7r „ t , > , . 1 \ 
1 2 > J n+v (2nx)J n -v\2nx) = —= 



and 



n=l 



v^F ^ n r(n + 1 + v)Y(n 

v n— 



2J2 Jn+ v [(2n + l)x]J n+ i- v [(2n + l)x] = -= ^ 

n=0 * *— " 



n! T(n 



+ 1 - z/) ' 
)(2x) 2n+1 



r(n + 1 + z/)r(n + 2 - 



(28) 



(29) 
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Figure 2. The Kapteyn series from equation I I29II for varying real v 6 [—10, 10] 
and for x = 1/2. 



which are illustrated in figures Q] and [21 respectively, for varying v and for x — 1/2. 
Numerically, the evaluation of infinite series is carried out as follows: First, a number 
of terms (usually 100) is summed directly; to a ccelerate the convergence of the sum, 



then for example Wynn's epsilon method (e.g., Brezinskl . l2000t lHammingl . Il986l ) can 
be used, which samples a number of additional terms (usually 100) in the sum and 
tries to fit a polynomial multiplied by a decaying exponential. Thus, the series are 
well approximated and the required computer time is kept moderate. 

Note that both equations (|28[) and (l29l) are valid for arbitrary v € C and for 
complex x € C with \x\ < 1/2. However, for > 1, both series attain extremely 
large values, depending on Sft(i/) and x. 

Nielsen then commented that for v = one has the particular cases 



1 + 2^ J n (2nxf = (1 - 4x 2 ) 



-1/2 



(30) 



and 



n=0 



J n [(2n + l)x]J n+1 [(2n ■ 



^ - Tx 



(1 - Ax 2 ) 



-1/2 



1 



(31) 



Note that the results shown here correct two misprints in Nielsen's results [his 
equations (30a) and (31)], which are: (i) Nielsen wrote J n - V [(2n + l)x] on the left- 
hand side of equation (j2"9"]l instead of J„ + i_^[(2n + l)x] and (ii) Nielsen included the 
term n = when summing the left-hand side of equation (|30p . 

In fact, however, a more general representation is possible if, instead of using 
Nielsen's series expansion of the Bessel function under the integral sign of equation (0, 
one were to write 



i r 

J2 n (2nxcos8) = — 

n Jo 



dip cos2n^cos(2nxcos , 0) 



(32) 
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then one sees immediately that the ability to evaluate a Kapteyn series is again reduced 
to the question of whether one can sum 

oo 

a n cos 2nip cos(2nx cos ip) 

71=1 

for particular choices of a n . Nielsen's choice of a n = 1 is just one example where the 
summation can be achieved. 

Perhaps of more general interest is to ask how reducible equation ([2]) is for 
arbitrary parameter values. Then, using equation (|27[) for each of the Bessel functions 
in K-i one has 



«* = E ««EE 



cn + b\ an+fj (-l) k (cn + b) 



21, 



2 2k k\T{an + /3 + k + l) 

/" + .'/) V""' i-i) r (/n + 3) 2 



n— — oo k— r— 

(33) 

2 2l >! r( 7 n + e + r + 1) ' v ; 

Unless fn+g = A(cn+6), where A is a constant, it is difficult to make further headway 
with equation ((331) - But when such is the case then one has 

K *= E '-EE 



n= — oo fc=0 r=0 



2 



(-l) fc (-l) r (cn + b^A 2 '- 

2 2fe 2 2r fc!H r(an + /3 + fc + l)r(7n + e + r + l)' 1 J 



where B n = a„A 7 " +e . Now, setting fc + r = to, one gets 



,x (a+7)n+e+/3 oo m / , 

i.= E^ EE 



f7l=0 /,■ = {) 

\m A 2m A — 2fc 



(-l) m A 2m A- 



(35) 

k\(m- k)\T(an + (3 + k + l)T(jn + e + l + m~ k)' y ' 

The representation of K2 in closed form (or at worst as an integral) then rests on the 
extent to which one can sum the various component sums occurring in equation (|35j) . 
One can write 

E o«. (36) 



where 



EE 



, \ 2m 

cn + 



fc=0m=0 

\m A 2m A —2k 



( _ 1} m A 2m A - ^ 



fc!(m - k)\ T(an + /3 + k + l)T(~fn + e + 1 + m - k) ' 

Note that the replacement of r by (to — k) causes Q n to be a single power series in 
{en + b)/2 which can be written as 

4 = Em (-l) m A 2m Rm (38) 
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R ™ = Y, 777 777 [ T ( an + (3 + k + l) T( T i + e + (m - k) + 1)] ~\ (39) 

£^ fc!(m — k)! 

The basic question is: under what conditions is R m expressible in closed form? If it 
is, then one can then determine the conditions under which Q n is expressible in closed 
form and so arrange values of B n so that K 2 is in closed form. It would seem that 
only for particular values of the parameters it is possible to effect closed-form results, 
as those for instance given by Nielsen [his equation (8)] 

f ) ^ = (/* + ,) r(i + M)r(i + „) £ (" + v + n l ) o* + * + 2n)-^ 

n=0 '' 



(I) 



x Jp+nMji + v + 2n)x] J u+n [(n + v + 2n)x] (40) 



3. Examples for K 2 series 

There are special cases of Kapteyn series of the second kind, which are often needed 
and which rely on other methods than those described above. One example is a 
Kapteyn series, which consists of J^(nz) for z € C with \z\ < 1, combined with a 
power of n in the form 

00 

K 2 = K 2 {n 2q ,1,0,1,0, «, 0,2,0) = ^n 2q Jl{nz). (41) 

71=1 

Two different distinctions can be made: (i) q < 0; (ii) q > 0; each for (a) integer 
q G IN; (b) arbitrary q G R. 



3.1. The case q < 
First, write p = — q so that 



K 2 



^ J n( UZ ) 



n 



2p 



(42) 



The general p rocedure is the following: use Bessel's equation J n (z), to show that (cf. 



Watsonl . 1 19661 equation 17.33), for consecutive indices p, the following two Kapteyn 
series of the first kind are related through the equation 

J2n(2nz) ; 1 . _> , \~» 

71=1 



z ■ 



n 



2p 



4(1 



J 2 „(2nz) 



n 2(p-l) 

n— 1 



(43) 



The second initial condition, i.e., the sum forp = 1 (cf. IWatsonl . [lQBfiL Sec. 17.23), 

E J 2 „(2nz) _ 
rT2 ~ T' 



(44) 



together with Meissel's <| 1892T ) investigation suggests one sho uld wri t e the 
Kapteyn series of the first kind as a polynomial in z 2k (see Watson . 19661 
Sec. 17.23). By evaluating the recurrence relation, equation (|4"3")l , it has been shown 
( Tautz and Dominici . 2010() that the Kapteyn series of the first kind can be expressed 
as 



E 



J 2n (2nz) 

n 2p 



yy fc y^ (-l) J+fc i 2(fc - y) 

k=i j 



(45) 




Figure 3. The functions ffc(p) for varying p with k £ {2, 3, 4, 5} as given through 
equation 1471 1 and, in explicit form, equation 1491 1. The vertical axis should merely 
illustrate the zeros of the functions §fc(p)i thus explaining why, for p integer, a 
finite power series is obtained in equation 1471 1, 



To obtain the corresponding Kapteyn series of the second kind, equation (|4"Tj) , 
one employs equation for /4 = u = n. Then one evaluat es equation (1431) with the 
argum ent 2nzcos(9 and integrates over 9, noting that (cf. iGradshtevn and Rvzhik , 
2000, Sec. 3.621) 



'^ C0S - 9 = M = M, „ eB , (46) 

(2n)\\ n\y/TT 

where (•)!! is the double factorial and where T(-) denotes the Gamma function. Hence, 
the result is 

KW-V,^V(i^ (47) 



BH t (48) 



with 

a_J 

oo, peR\B 

and is valid for arbitrary p £ R < 0. 

Equation (|47[) has some important features. For integer p, a sum with p{p+ l)/2 
terms is obtained, whereas, for non-integer p, an infinite power series occurs. But even 
in that case, it is advantageous to exchange one infinite series (the Kapteyn series) by 
another infinite series (a power series), because the convergence behavior of a power 
series is better understood, thus allowing for a more reliable estimate of the number 
of terms needed to obtain a desired accuracy. 

The reason for the distinction between a finite/infinite power series is that, when 
expanding the coefficients of the powers z 2k , one finds non-algebraic functions £/. (p) 
of the form 

UP)= -^ + 4~ (1+p) 
16 
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e 3 (p) = 7 gg(5-2 7 - 2 P + 3 5 ~ 2 f) 



Up) 
Up) 



i 



18 432 
1 



(-7 + 2 13 ~ 4p + 7 • 2 7 ~ 2p - 9 4 ~ p ) 

(42 - 2 21 ~ 4p - 3 ■ 2 13 ~ 2p + 5 9 - 2p + 9 6 ~ p ) , (49) 



2 949120 

each of which has zeros at the first integers, i. e., £fc = for p = 1, . . . , k — 1. The 
functions £fc(p) are illustrated in figure |3]for varying p and for k G {2, 3, 4, 5}. 

3.2. The case q>0 

Even though structurally similar, the opposite case with q > results in a completely 
different b ehavi or when evaluating the Kapteyn series in terms of a power series. In 
Dominici ( 201 lh . Kapteyn series of the form 



K 2 (z, q) = K 2 (n 2q , 1, 0, 1, 0, 2z, 0, 2z, 0) = J 2 (2nz) , (50) 

71 = 1 

were investigated and it was found that 



(51) 



*.c«)-tfwt(- i )*cr)("-*) w 

n=l Lv'v fe = o v 7 

for g ^ and |z| < 1/2. The result in equation (|5"Tj) allows one to compute K 2 (z,q) 
numerically for any q ^ 0,; however, sometimes it is difficult to use it to get closed-form 
expression s. Here, a differ ent approach will be introduced using differential operators. 

From IWatson ( 19661 equation 5.4 (4)), it is known that the function y n (i) = 
J 2 (e*) satisfies the differential equation 



d 3 y n 



•4(e J 



di» > -\r-n-)-fr + 4B y n = 0. 



By changing variables to 2nz — e* in equation (|52|) , one obtains 



d 3 



z~-r^r + %z—z + (l6n~z 
dz- 3 dz z 



2 - 2 4n 2 + l) — + 16n 2 z 



J 2 (2nz) = 0, 



or 



dz 3 dz 2 dz 
Introducing the function 

g q (z, n) — n 2q J 2 (2nz) 

one has 

2 d 3 d 2 d , 



4(1 -4z z ) — - 16 



dz 



n 2 J 2 (2nz) . 



dz 2 



dz 



4(1 -4z 2 ) — - 16 
dz 



for all n £ IN*, i.e., all positive integers. 

Thus, it follows that the function K 2 (z,q) satisfies 



dz 3 dz 2 dz 



K 2 (z,q) 



4(1 -4z 2 ) 16 

dz 



K 2 (z,q + 1), 



while equation (|2|) gives for the initial condition 

1 1 



K 2 (z,0) 



2 2Vl-4z 2 ' 



(52) 

(53) 
(54) 
(55) 
(56) 

(57) 
(58) 
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K 2 (z,q) 

10 7 r 




0.02 0.05 0.10 0.20 0.50 

Figure 4. The first functions K2(z, q) for varying z S [0, 1/2] with q £ {1, 2, 3, 4} 
as given through equation H62I I and, in explicit form, equation H61I I. The order q 
varies from right to left, i.e., the solid line shows K2(l,z) while the dotted line 
shows if 2 (4, z). 



Since J n (0) = for all n = 1, 2, ... , one has ^2(0, q) = 0. Solving equation (|57|) for 
i4T 2 (z, <? 4- 1), one obtains 



K 2 (z,q+1) = 



1 



dw 



Vl-4^ 2 7 n/1 - 4w 2 V dw3 



where one must be careful that \z\ < 1/2 to ensure that the integral is convergent. 

Integrating equation (|59l) by parts, one obtains the final result for the series K 2 
in the form of a the recurrence relation, which reads 



K 2 (z,q+1) 



(l-8z 2 ) d 2z 2 (l + 2z 2 



4(l-4z 2 )dz 2 4(l-4z 2 ) 2 dz (1 - 4z 2 ) J 
4 



K 2 (z,q) 



u; (1 + 10w 2 + 4w 4 ) 
du; — i if 2 (w, ?). 



(1 - Aw 2 ) 



7/2 



(60) 



Using equation (|60|) and (|58|) . it is straightforward (albeit tedious) to compute the 
explicit expressions for the first orders of the series K 2 , yielding the relations 

^2 



K 2 (z,l) = 
K 2 (z,2) = 
K 2 (z,3) = - 



(1 + ^ 2 ) 
(l-4z 2 ) 7/2 

z 2 (l + 37z 2 + 118z 4 + 27z 6 ) 



(1 -4z 2 ) 13/2 

(1 + 217z 2 + 5036z 4 + 23 630z 6 + 22 910z 8 + 2250z 10 ) 



(1 - 4z : 



,19/2 



K 2 {z, 4) = z 2 (385 875z 14 + 7 119 756z 12 + 15 359 862z 10 + 8 635 5782 



1 515 705z 6 + 80130z 4 + 973z 2 + l) (l - 4 



2 x-25/2 



(61) 
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which are illustrated in figure 2J 
In general, one has 

K ^ q)= { i-2^ ' N 4 (62) 

where P q (z) is a polynomial of degree 2q— 1. The structure of the polynomials -P<j(z) 
is quite complicated and will be analyzed in a forthcoming paper. 

4. Discussion and Conclusion 

In this article, the general features and characteristics of Bessel function series were 
investigated. Special emphasis was focused on Kapteyn series, which appear in 
many applications of theoretical physics and mathematics, such as radiation and 
optimization problems. In their original form, with the index of summation appearing 
in both the index and the argument of the Bessel function(s) involved, the convergence 
of such series is, in general, unclear. Therefore, it is appropriate and necessary to 
undertake every effort of rewriting such sums in terms of, at worst, infinite power series 
or (double) integrals, the convergence of which can be estimated more reliably. More 
importantly, in many cases it has proven possible to find closed analytical expressions 
for Kapteyn series of both the first and the second kind. This is indispensable for cases 
where Kapteyn series constitute only part of large mathematical expressions that have 
to be dealt with numerically. 

However, the quest for closed-form expressions of Kapteyn series has mostly been 
limited to special cases that have arisen in specific problems such as those listed in 
Sec. [1] Often one notes that a similar procedure proves useful for different forms of 
Kapteyn series, even in such cases where the summation coefficients are considerably 
diverse. But no general answer has been found to date to the following problem: For 
which parameter regimes of the coefficients do the Kapteyn series have a closed-form 
expression? 

It was the aim of the present article to shed some light on that question. 
Starting from the most general form of Kapteyn series of the second kind, i.e., 
involving the product of two Bessel functions, the Kapteyn series was decomposed 
into series over trigonometric functions or, more generally, algebraic expressions 
involving gamma functions. The ability to sum such series depends on the precise 
choice of the parameters. However, it has been shown that the likelihood for any 
analytical tractability is increased if both Bessel functions are of equal order, i.e., 
(an + ft) — (771 + e) in equation ©. 

Tw o specific examples with applications, for instance , in cosmic ray diffusion 
theories ( Tautz and Dominici , 20101 Shalchi and Schlickeiseii 12004 iTautz et al.ll2006l) 



were illustrated, where the summation coefficients are simply powers of the summation 
index. It has been shown that, depending very sensitively on the parameters chosen 
for the summation coefficients, power series with a finite/infinite number of terms are 
obtained, where the relative magnitude of each term can now easily be estimated. 

Future work should, presumably, concentrate on the application of one or more 
of the above methods (or indeed combinations of the methods) to determine to the 
maximum extent possible the broadest range of conditions for summability of the 
general Kapteyn series. Of course, new methods for attempting such summations are 
welcome and it would be highly interesting to see any such ideas that would add to 
the capability to effect Kapteyn series summations. Another interesting question is to 



Summing General Kapteyn Series 



14 



decide if for every Kapteyn series there exists an annihilates differential operator that 
can be used to obtain structural relations between series of different orders. 
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